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(57) Abstract 



The present invention relates 
to the use of time gated scattered 
light, for determining the location 
and composition of material within 
various organs of the human body. 
The systems and methods of the 
present invention provide for medi- 
cal imaging in three dimensions of 
internal body structures for diagnos- 
tic purposes. Stationary and scan- 
ning fiber optic systems can be used 
deliver laser radiation onto the ob- 
ject to be imaged and to collect the 
gated information for image re- 
construction. 
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OPTICAL IMAGING USI NG TIME GATED SCATTERED LIGHT 

Related applications 

This is a Continuation-in-Part application of U.S. 
Application Serial No. 08/549,623 filed on October 27, 
5 1995, which is a continuation-in-part application of U.S. 
Application Serial No. 08/391,209, which was filed on 
February 21, 1995, the contents of these applications being 
incorporated herein by reference in their entirety. 

Background of the Invention 

10 Th e present invention relates to medical and diagnos- 

tic imaging systems and procedures. Clinical procedures 
currently employ a number of systems for locating, imaging 
and diagnosing various structures within the human body. 
These include x-ray computer tomography, ultrasound, and 

15 magnetic resonance imaging, among others. These systems 
are used to detect morphologic abnormalities associated 
with specific diseases or conditions in various body or- 
gans . 

In the case of x-ray computer tomography, for example, 
20 a number of projection data are taken sequentially at 

different angles and the data are used to reconstruct an 
image of the object being scanned i.n three dimensions. 
Thus, an x-ray tomography system solves an inverse problem 
for the x-ray opacity of body tissues using measurements of 
25 the amount of radiation absorbed from many beams transmit- 
ted at a variety cf angles. This procedure is based on a 
number of assumptions including that the intensity of the 
x-rays diminishes across the distance traversed at a rate 
proportional to the intensity of the beam, that the absorp- 
30 tion coefficient depends on the type of tissue along the 
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various beam paths and that this non- linear problem can be 
solved based upon certain approximations including that a 
linear set of equations is an accurate representation of 
the problem. Of critical importance to x-ray tomography is 
5 that corrections for scattering are relatively simple. 

Others have sought to use the "diffusion approxima- 
tion" to represent the scattering of optical radiation for 
medical imaging applications. The diffusion approximation 
involves the detection of incoherent photons and the analy- 

10 sis of the resulting spectrum. The main problem with this 
approach is poor spatial resolution which limits the use- 
fulness of this method in medical imaging applications. 
Others have sought to use so-called "ballistic" photons 
which travel the shortest path through the medium and are, 

15 for the purposes of this application, "non-scattered" 
photons . 

A continuing need exists, however, for further im- 
provements in the field of tomographic imaging for medical 
applications including enhanced resolution, reduced cost 
20 and complexity, and improved diagnostic capability. 

Summary of the Invention 

The present invention relates to the use of time 
resolved elastic and inelastically scattered light for 
locating, imaging and diagnosing structures within organs 

2 5 of the human body and devices used for the delivery and 

collection of opcical radiation for this purpose. In 
particular time-resolved photon migration is used for 
medical imaging, including optical methods for localizing 
lesions within the body. Since biological tissue is highly 

3 0 scattering, the problem is one of imaging an object embed- 

ded in a turbid medium. Most existing techniques use 
differences in the absorption or elastic scattering proper- 
ties between the embedded object and its surroundings. In 
many cases of medical interest, however, the resulting 
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contrast is relatively small, which severely limits the 
sensitivity and resolution. 

Fluorescence spectroscopy studies of human tissue 
indicate that a variety of lesions show distinct fluores- 
5 cence spectra compared to those of normal tissue. Thus, 
intrinsic tissue fluorescence can provide enhanced con- 
trast, as well as diagnostic-histochemical information. In 
addition, exogenous dyes, many of which are known to fluo- 
resce with high quantum yield, have been shown to exhibit 

10 selective uptake in neoplastic lesions. Use of such agents 
provides fluorescent markers with high quantum yields, and 
are used to locate embedded lesions in the breast, brain 
and other organs . 

In a preferred embodiment of the present invention 

15 time-resolved optical tomography, and fluorescence and 

Raman spectroscopy, can be used separately or in combina- 
tion to provide both spatial and chemical information about 
embedded objects in tissue. By measuring and analyzing the 
early portion of the fluorescence signal from embedded 

20 lesions for example, which rises rapidly and is not sensi- 
tive to fluorescence lifetime, precise timing information 
and hence accurate spatial resolution of embedded lesions 
can be obtained. The rising edge of the fluorescence 
signal is generally over a period of 100-200 picoseconds or 

25 less. The rising edge of the Raman signal is less than 100 
picoseconds. A streak camera can be used as a multichannel 
time-resolved detector to measure both the rising and 
decaying periods of the spectrum and to obtain images of 
embedded fluorescent objects in a single measurement or in 

30 a sequence of measurements at one or more angles relative 
to the tissue. 

More specifically, a preferred embodiment collects 
scattered radiation in the interval between 0 and 1500 
picoseconds, and preferably in the range of 0 to 500 pico- 

35 seconds after irradiation of the tissue, and based upon a 
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comparison of this data with a non-diffusion representation 
of light which has been scattered by the tissue. This non- 
diffusion representation emphasizes the "almost straight" 
trajectories of early arriving photons to provide images of 
5 internal bodily structures with improved spatial resolution 
in the range of 1-3 mm or less. This representation can 
include a diffusion component for later arriving photons in 
each collection period which exhibit characteristics more 
accurately represented by both diffusion and non-diffusion 

10 characteristics. 

Time resolved Raman scattering of tissue involves the 
detection of early arriving photons arising from molecular 
vibration in the objects being measured. The vibrational 
bands can be assigned to individual molecular groups so 

15 that information about molecular content as well as posi- 
tion can be provided. 

Optical collectors positioned about the object to be 
imaged are used to collect fluorescence and/or Raman light 
scattered by the tissue in response to laser irradiation. 

20 The fluorescence and Raman data can be used to locate and 
image embedded lesions in three dimensions as well as 
provide information regarding the chemical composition of 
such lesions. This information can be used to identify 
such lesions as normal or abnormal, cancerous or precancer- 

25 ous etc. Further details regarding the use of induced 

fluorescence of tissue to diagnose cancer or precancerous 
lesions can be found in U.S. Serial No. 08/219,240 filed on 
March 29, 1994, now U.S. Patent No. 5,421,337, the entire 
contents of which is incorporated herein by reference. The 

30 use of laser induced Raman spectroscopy of tissue for 

diagnosing various diseases and conditions is described in 
greater detail in U.S. Serial No. 08/107,854 filed on 
August 26, 1993 and having an international filing date of 
January 17, 1992, and entitled Systems and Methods of 

3 5 Molecular Spectroscopy To Provide For The Diagnosis of 
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Tissue, the entire contents of which is incorporated herein 
by reference. 

Another preferred embodiment utilizes back-scattered 
fluorescent and/or Raman data to provide measurements of 
5 depth and lateral position of lesions in tissue. Fiber 
optic probes used for back scattering measurements are 
described in conjunction with the above referenced incorpo- 
rated applications. The time resolved methods can be used 
or, alternatively, frequency domain analysis of the ac- 

10 quired spectrum can also be employed. In this embodiment, 
it is desirable to modulate the laser above 100 MHZ, and 
preferably between 500 MHZ and 1000 MHZ to provide the 
temporal resolution for accurately imaging embedded ob- 
jects. A computer is programmed to transform the data 

15 using known transformation techniques and the data can be 
represented in three dimensions with amplitude plotted as a 
function of frequency and time. 

A preferred method of using the invention to diagnose 
tissue involves the insertion of a fiber optic probe or 

20 endoscope into a body lumen, collecting data, determining 
the location of material to be biopsied, and performing a 
biopsy on the identified material. The proximal end of the 
device can be coupled to a streak camera or other suitable 
detectors such as a change coupled device. Additionally, 

25 biopsied samples can also be measured and analyzed using 
the methods described herein. Alternatively, an optical 
biopsy can be performed in vivo without the need for tissue 
removal using the methods set forth herein. 

Other preferred embodiments of the invention relate to 

30 the imaging of subsurface lesions using scanning or sta- 
tionary optical delivery and a collection systems. The 
particular system selected will depend upon the particular 
type of lesion or tumor, the depth of the target tissue, 
the resolution desired, and the time necessary to collect a 

35 sufficient number of photons from different views of the 
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target sufficient to construct the desired representation. 
In a first embodiment a single ring of collection fibers is 
scanned across a target. In a second embodiment, a partial 
ring is used. In a third embodiment a number of rings are 
5 stacked to provide three dimensional information. in a 
fourth embodiment, a collection system is positioned rela- 
tive to the target which collects over a relatively small 
field of view while the source of radiation scans across 
the target either mechanically or optically. 

10 The above and other features of the invention includ- 

ing various novel details of construction and combinations 
of parts will now be more particularly described with 
reference to the accompanying drawings and pointed out in 
the claims. It will be understood that the particular 

15 system and methods embodying the invention are shown by way 
of illustration only and not as a limitation of the inven- 
tion. The principles and features of this invention may be 
employed in varied and numerous embodiment without depart- 
ing from the scope of the invention. 

20 Brief Description of the Drawing s 

Figure 1A is a schematic diagram of a tomographic 
apparatus where 0 is an embedded fluorescent object; D is a 
fast photodiode; F are collection fibers; HF is a holo- 
graphic notch filter. 
25 Figure IB schematically illustrates a circular 

multifiber collection geometry. 

Figure 2 is an 8-channel time - resolved fluorescence 
signal observed by the streak camera where the line con- 
nects t 1/2 points of the eight channels. 

Figure 3A is a graphical representation where fluores- 
cence signal detected by each individual fiber surrounding 
the object containing a single embedded object. Fits using 
a triangulation procedure are plotted as dashed curves. 
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Figure 3B is graphical representation similar to that 
of Figure 3A with a lower mean free path. 

Figure 4 is a graphical representation where t 1/2 is 
the fluorescence signal detected by each individual fiber 
5 surrounding a medium containing two embedded objects. 

Figure 5A illustrates a schematic diagram showing 
instrument components where CD: cavity dumper; MCP/PMT: 
microchannel plate/photomultiplier tube; CFTD : constant 
fraction time discriminator; TAC: time-to-amplitude con- 
0 verter; and PC: personal computer. Solid lines represent 
optic signal paths and dashed electronic signal paths. The 
insert shows the arrangement of excitation and collection 
probes, and the sample geometry. 

Figure 5B schematically illustrates an endoscopic 
5 backscattering geometry for a fiber optic probe. 

Figure 6A graphically illustrated a fluorescence 
emission from rose bengal (RB) dye cell at four different 
depths . 

Figure 6B is a graphical illustration similar to that 
0 of Figure 6A with an rhodamine (R6G) dye cell. 

Figure 7 graphically the depths of a fluorephore vs. 
t 1/2 for two different scattering coefficients. 

Figure 8 graphically illustrates lateral location of a 
fluorophore using ultra fast time-gating for two different 
5 scattering coefficients. 

Figure 9 graphically illustrates the probing two 
identical fluorophores at different depths. 

Figures 10A and 10B graphically illustrate two differ- 
ent fluorophores detected at 620 nm and 670 nm respective- 
0 ly . 

Figure 11 graphically illustrates results with Raman 
scattering from /3-carotene cell at three depths. 

Figure 12 graphically illustrates time resolved mea- 
surements of a turbid medium. 
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Figures 13A and 13B graphically depict the half of the 
maximum width of path distribution r(a), and relative 
intensity of emerging light Kb) as a function of time for 
3 cm of turbid medium representing breast tissue (fi' s = io 
5 cm" 1 ) found from path integral formulation of transport 
equation. 

Figure 14 illustrates detection of the screen edge in 
a 5.5 cm slab of turbid medium in' s = 0.45 cm" 1 ) assuming a 
uniform distribution of the trajectories. Lateral position 
0 of the edge is d=l cm. Theory (solid lines) vs. experiment 
(dots) , where different curves represent time windows at 
t=310 ps, 450 ps, 500 ps, and 600 ps . 

Figures 15A-15E schematically illustrate preferred 
embodiments for the delivery and collection of radiation 
5 for imaging of tissue. 

Detailed Description of the Invention 

A schematic diagram of apparatus in accordance with 
the invention is presented in Figure 1A. The system can 
use -150 fs excitation pulse generated by a Coherent Mira 

0 900 mode-locked Ti : sapphire laser 10 pumped by a Coherent 
Innova 400 multiline argon ion laser, and a streak camera 
detection system 22 consisting of a temporal dispenser 
C1587, synchroscan streak unit M1955 and tuning unit M1954 
available from Hamamatsu Photonic Systems Corporation. A 

5 photomultiplier detector can also be used. The wavelength 
in the present embodiment is 750 nm, the repetition rate 76 
MHz, and the average power 1.5 W although other wavelengths 
rates and power levels can be used. A small portion of the 
excitation light, deflected by a quartz plate to a fast 

0 photodiode 20 (D) , was used as the optical triggering 

signal. Eight 100 ^m cere diameter optical fibers (F) were 
used to collect the fluorescence light and transmit it to 
the entrance slit of the streak camera using a fiber optic 
cable 18. The excitation light can be delivered to the 
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sample along fiber optic cable 14. An endoscopic probe 
used in conjunction with the invention can incorporate both 
delivery and collection fibers. A Kaiser 752 nm holograph- 
ic notch filter 26 (HF) and two 780 nm long-pass filters 
5 were used in front of the streak camera to completely 

remove the excitation light. The system resolution, 10 ps, 
was determined by the intrinsic response of the streak 
camera, and temporal dispersion through the optical fibers, 
and optical trigger jitter. 
0 Signals can be collected by 8 fibers evenly positioned 

around the circumference of the medium in the plane con- 
taining the fluorescence object (0) as shown in Figure IB. 
Other numbers of fibers and fiber geometries can also be 
used. Differences in length among the fibers can cause 
5 uncertainties in the zeroes of time, but these can be 
calibrated out. In this arrangement, as shown is Figure 
IB, the fluorescence return signal to be received earliest 
by the distal end of fiber 2, which is closest to the 
fluorescent object, and latest by fiber 6, which is fur- 
thest away. The proximal ends of these fibers can be 
arranged in a line and imaged onto the streak camera slit. 
The streak camera was used as a both time-and spatial- 
resolved multichannel detector. Signals from different 
fibers appeared at different vertical locations on the CCD 
array of the streak camera, and the temporal information 
was displayed horizontally. The detected data can then be 
processed by computer 24 to provide the desired time-depen- 
dent representation and image of the tissue or features 
therein and identify the fluorescent or molecular compo- 
nents of interest. 

A typical 8-channel time -resolved fluorescence signal 
is shown in Figure 2. The rising edge of each curve con- 
tains timing information which can be used to extract 
spatial information. The exact timing of these rising 
edges can be determined in several ways. In the present 
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method the time at which the signal reaches half -maximum of 
that channel, t 1/2 was selected at the representative time 
for each curve. In Figures 3A and 3B we plot t 1/2 for each 
collection fiber for a single imbedded object for the cases 
5 of 1.8 and 0 . 9 mm scattering mean free path (mfp, l//i„) , 
respectively. Note that ^ s and /i a are the scattering and 
absorption coefficients of the medium. The collection time 
is 5 minutes in the example, but can be longer or shorter 
depending on the object to be imaged. The local minimum 

10 indicates the position of the object. Similar data were 
obtained for mfp's between 0.6 and 3.6 mm. (Note that the 
diameter of the sample is equivalent to 20-120 mfps) . 

In order to extract the spatial location of the embed- 
ded fluorophore, in this procedure it is assumed the time 

15 delay for the early photons to be proportional to the 

distance traveled. In our cylindrical geometry, this can 
be expressed as 

r n °yi? 2 +r 2 -2ri?cos(e-e n ) , 

with R the radius of the sample, 0 n the angular location of 
the nth channel, and (r, 6) the polar coordinates of the 

20 embedded fluorophore. The spatial location of embedded 
objects can be obtained by fitting the experimental data 
(e.g. Figures 3A and 3B) using Eq. (1). Typical fits are 
shown as dashed lines in Figures 3A and 3B. 

To illustrate the measurement of two embedded objects 

25 in which the scattering mfp was 1.8 mm, Figure 4 plots t 1/2 
for each collection fiber with two objects present. Each 
local minimum represents a single fluorescing object. Use 
of more than eight fibers improves the resolution particu- 
larly when measuring a number of fluorescing objects. The 

30 inverse problem in this case is more complicated and re- 
quires an accurate theoretical model of the photon migra- 
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tion process in a turbid medium described in greater detail 
below. 

This geometry is applicable to many human organs, such 
as brain and breast. Furthermore, this method can be 
5 extended to other geometries, as well. Note that with 
fibers positioned in different planes, objects can be 
localized in three dimensions. In addition, by observing 
fluorescence at multiple wavelengths, one can obtain histo- 
chemical information. 

10 The streak camera photocathode -window (s-20/UV, spec- 

tral response from 200-850 nm) can be more nearly matched 
to the fluorescent wavelength or wavelengths depending on 
the fluorescence wavelengths of interest. In addition, the 
system throughput can be improved by replacing the above 

15 fibers with single-mode optical fiber bundles of larger 

diameter, e.g., 1 mm, without loss of temporal resolution. 
Improving these factors increases the signal-to-noise ratio 
by a factor of several hundreds, which permits reduction of 
the excitation power and/or collection time. The localiza- 

20 tion accuracy comparable to that of the data presented here 
can be obtained from a 5 mm diameter fluorescent object 
containing hematoporphryn derivative (HpD) at clinical 
concentrations using excitation power of a few tens of mW 
and less than one minute accumulation time. 

25 In another preferred embodiment the rise- time of 

fluorescence provides precise timing information, hence 
accurate spatial resolution, of an embedded object in 
tissue. Electronic excitation occurs on the optical time 
scale; and the rise time, which is not affected by fluores- 

30 cence decay, is on the order of 10" 13 seconds. Raman scat- 
tering, another potential source of high contrast, exhibits 
comparably rapid rise times and can be used for precise 
timing. Note that these spectroscopic techniques carry the 
additional potential advantage of providing histochemical 

35 information about the object, to be imaged. Thus fluores- 
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cence emission and Raman scattering are used to obtain 
time-resolved and three-dimensional spectroscopic informa- 
tion of lesions through overlying tissue. 

A preferred embodiment can use a few-ps pulses for 
5 excitation and time-resolved single-photon counting for 
detection. A wavelength suitable for this application be 
generated using a Nd : YAG laser 30 to provide radiation at 
570 nm, the repetition rate is 1 MHz and the average power 
is 3 0 mW. The apparatus is illustrated schematically in 

10 Figure 5A. A 10 nm bandpass filter (centered at 610 nm) 

and a 600 nm long pass filter in front of the detector com- 
pletely remove the excitation light. The temporal resolu- 
tion of the system, about 80 ps , is mainly due to the 
transit time spread of the PMT. This is adequate to re- 

15 solve the photon migration signals in our photon samples. 
Tube 3 8 is used with a microchannel plate. A constant 
fraction time discriminator 44 is coupled to tube 38 and 
the time to amplitude converter 46. An endoscope 36 can be 
used for delivery and/or collection in backscattering probe 

20 procedure for in vivo applications {figure 5B) . A diode 4 0 
and second CFTD 42 can be used to improve collection prior 
to analysis on computer 48. 

In our single-ended probing geometry, a collection 
fiber is positioned adjacent to the incident beam on the 

25 same surface of the tissue, and the fluorescence material 
is located at a fixed distance within the tissue. In this 
geometry the earliest arriving fluorescence photons emitted 
from the fluorescing material or lesions returns to the 
collection fiber sooner than those from a second lesion B, 

30 located at a greater depth, or a third lesion C, laterally 
displaced from A. This illustrates the use of a fiberoptic 
probe in determining the location of the fluorescent object 
in three dimensions by observing the time-resolved fluores- 
cence signal. Ultrafast time-gating can further improve 

35 spatial resolution. 
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Fluorescence can be used to provide time-of -flight 
signals for ranging, despite the long fluorescence life- 
times,, typically on the order of a few ns, relative to the 
photon migration time of interest. To illustrate this, we 
5 compared the time-resolved signals emitted from dye cells 
embedded in a polystyrene bead suspension containing either 
RB (in water lifetime 130 ps) or R6G (lifetime 3.9 ns) . 
The dye cell depth was changed from 3 to 6 cm in 1 cm 
increments. The scattering mean free path (mfp) was 4 mm. 

10 The results are plotted in Figures 6A and 6B . In spite of 
the fact that the fluorescence lifetimes of RB and R6G 
differ by a factor of 30, the rising edges of these two 
sets of curves are extremely similar. In fact, they close- 
ly represent the rising edge of the photon migration 

15 curves. 

The depths of a fluorescent object can be easily 
probed in the backscattering geometry by looking at the 
earliest-arriving fluorescence photons, even under strong 
scattering conditions. The time-of -flight for these earii- 

2 0 est arriving photons is proportional to twice the depth of 
the fluorescent object. As expected, fluorescence emitted 
from a deeper object arrives at the detector later. In 
addition, increased scattering also delays the arrival 
time, as shown in the plot of depth vs. time, of Figure 7. 

25 The time at which the signal reaches half maximum, t 1/2 is 
used as the representative time for each curve. The spa- 
tial resolution is determined by both the intrinsic factor, 
due to the statistics of photon migration, and extrinsic 
factors such as the sensitivity of the detector and signal - 

30 to-noise (S/n) level. With five-minute accumulation time 
the current set-up can resolve a depth of 1 mm, even at a 
distance of 100 mfp's. 

To obtain lateral information about an embedded ob- 
ject, we measured the time-resolved fluorescence signals as 

35 a function of the lateral position of the dye cell at depth 
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of 5 cm. Based on simple geometrical considerations, we 
expect lateral resolution to be poorer than depth resolu- 
tion. However, lateral resolution can be improved by using 
an ultrafast time window for the earliest arriving photons, 
5 which is similar to the typical time-gate transillumination 
imaging experiments. The optimum time gate was chosen 
based on considerations of both resolution and S/N ratio. 
Typically, we used a 24 ps time window. Figure 8 plots the 
time gated intensity vs. the lateral displacement of the 

10 fluorescence dye cell at 25 and 100 mfp's, respectively. A 
lateral resolution of 1 cm is obtained with this time win- 
dow, however, much greater resolutions can be obtained. 

Measurements were made with two embedded cells con- 
taining the same fluorescence dye. The two objects were 

15 positioned 1 . 8 cm apart at depths of 4 . 5 and 5 cm, respec- 
tively. Using lateral displacement and time delay as 
variables, we can construct a contour map of the time- 
derivative of the time-resolved intensity (Figure 9) . The 
time-derivative, a measure of how fast the time-resolved 

2 0 signal evolves, reaches a maximum at approximately t 1/2 , 

which in turn provides the depth information. Although the 
signal from the shallower object dominated, the 3-d posi- 
tions of both objects can be ascertained. This is because 
we can use different time-gating to provide different depth 

2 5 measurements inside the tissue medium. Note that although 
the lateral dimension represents the actual image of the 
object, longitudinally (i.e., temporally) it only provides 
information about the top surface of the fluorescent mate- 
rial. Also note that the best lateral resolution is ob- 

30 tained at the earliest time gate, i.e., where the contours 
begin along the time axis. 

If the two embedded objects have distinct spectroscop- 
ic features, imaging capability can be enhanced. This was 
demonstrated by measuring the fluorescence intensity from 

35 two cells (similar geometry as in Figure 9) containing S640 



WO 96/26431 



PCT/US96/02394 



-15- 

and HIDCI, respectively, at 620 and 670 nm (Figures 10A and 
10B respectively) by placing a spectrometer in front of the 
PMT. Compared to Figure 9, signal interference between the 
two objects is reduced, and more importantly, the chemical 

5 identities of the objects can be obtained. 

A Raman scattering cell containing /3-carotene having a 
Raman vibration at 1157 cm" 1 was studied. To establish 
that the detected signals were from Raman scattering rather 
than fluorescence from either ^-carotene or sample impuri- 
10 ties, in a measurement without the scattering medium, the 
time-resolved Raman scattering was identical in shape to 
that of the laser light (determined by instrumental resolu- 
tion) , obtained by deflecting part of the laser beam into 
the fiber probe and removing the filters. This is consis- 

5 tent with the instantaneous nature of Raman scattering. 
The possibility of laser light leakage was excluded by 
replacing the Raman cell with a Raleigh scattering cell 
containing KI powder, which does not have Raman bands above 
300 cm" 1 . As a further check, the emission spectrum from 

0 the Raman sample was collected using a standard fluorimeter 
with the same 570 nm excitation. Distinct Raman peaks at 
605, 610 and 624 nm, which correspond to the known Raman 
shifts of 1008, 1157 and 1516 cm" 1 for /^-carotene, were 
observed with minimum background. Figure 11 displays the 

5 time-resolved Raman signal for cell depths of 2 , 3 and 4 cm 
in a scattering medium with 7 mm mf p . 

Thus time-resolved fluorescence and Raman emissions 
can be used as a probe to provide accurate information 
about the position of an object such as a lesion embedded 

0 in a turbid medium such as tissue. Time resolved fluores- 
cence can be used effectively despite the long fluorescence 
decay time. The time evolution of the signal is determined 
by both the time-course of the photons migrating through 
the scattering media and the decaying fluorescence. These 

5 measurements establish that the early portion of the signal 
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rises rapidly and is not sensitive to the relatively long 
fluorescence lifetime. The early portions of both the 
fluorescence and Raman signals thus represent the actual 
time-of -flight of photons traversing the scattering medium. 
5 In the case of fluorescence, this is due to the fact that 
electronic excitation from the ground state to the excited 
states is extremely rapid, on the order of 10" 13 seconds, so 
that the rising edge of the fluorescence lifetime curve is 
effectively a step function. The resulting signal, a 

10 convolution of the distribution of the time-of - flight of 
photons migrating to the detector and the fluorescence 
decay curve, is mainly determined by that of the time- of - 
f light curve. In addition, these earliest arriving photons 
undergo minimal scattering, resulting in high spatial 

15 resolution. In other words, by measuring and analyzing the 
earliest-arriving photons, this serves to minimize the 
uncertainties due to both diffusive scattering and the 
finite fluorescence lifetime. 

However, unlike many existing trans- illumination 

20 methods, a backscattering geometry with single-ended detec- 
tion is employed. In this configuration, the isotropic 
fluorescence emission or Raman scattering events serve to 
reverse the direction of nearly straight photon paths. The 
earliest signals observed in this case are proportional to 

25 twice the distance from the sample surface to the location 
of the embedded objects, which in our experiments are 
composed of molecules having unique fluorescence or Raman 
scattering characteristics. Thus, the arrival time of the 
earliest photons can provide the depth information of the 

30 embedded objects, as well as the usual 2-D localization. 

In addition, lateral spatial resolution can be optimized by 
using a short duration time-gate to select the earliest 
arriving photons. Also, note that in a turbid medium, a 
backscattering geometry cannot be used with a technique 

35 based on elastic scattering or absorption, because photons 
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backscattered by the medium will mask the signals arising 
from the embedded object. 

With the technique presented in this paper, either 
native tissue chromophores or exogenous dyes can be used to 
5 detect embedded lesions. In the former case, both biochem- 
ical and spatial information about the lesion are provided 
simultaneously. The contrast between diseased and normal 
tissue can be enhanced by utilizing these properties in- 
stead of absorption or elastic scattering, and diagnostic 

10 information can be obtained. Exogenous dyes, which are 

known to exhibit fluorescence with high quantum yield, can 
also be used as a source of contrast. For example, a 
number of known methods use selective uptake of photosensi- 
tizing agents, such as hematoporphyrin derivative, in a 

15 neoplastic lesions. Use of such agents provides fluores- 
cent markers with high quantum yields which can serve to 
locate embedded lesions in the breast, brain or perhaps 
other organs. The dyes preferably have an excitation 
wavelength in the range of 600-630 nm and fluorescence in 

20 the range of 680-720 nm. 

Multichannel detection can be used to extend the 
capability of the technique presented here in several ways. 
An optical fiber array can be utilized to obtain 3-D infor- 
mation in a single measurement by displaying physical 

25 position along one axis and time along the other. Further 
more, with the aid of a spectrograph, real-time spectrosco- 
py and optical tomography can be achieved by simultaneously 
displaying spectral and temporal information. 

To understand the differences between prior methods 

30 and the present invention it is helpful to consider a 

system in which photons propagating through a turbid mediurr. 
from point A to point B can be described using statistical 
considerations. A photon originating at point A is elasti- 
cally scattered multiple times, each time being deflected 

35 into a particular angle with a well-defined probability, 
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thus forming a trajectory within the medium. Each trajec- 
tory has a particular probability. By calculating this 
probability and then summing over all possible trajecto- 
ries, one obtains the probability for the photon to travel 
5 between two points in the medium. 

The probability for a photon to traverse the scatter- 
ing medium from point A with radius-vector r A to point B 
with radius-vector r B in a time interval T can be written 
as a path integral : 

T 

P(r A ,r B ,T) =jl>r(t)exp{- ^ 1 — ft ( t) 2 dt)j[f ( t) ] , (2) 

T 

where the function J[r(t) ] =JdQ ( t) exptiji.JdtQ ( c) (f (t) 2 -l)} 

o 

10 insures that photons propagate at the speed of light in the 
medium (c=l) at every point along the trajectory, and 

N 

Dr(t)=Y[dr n - O) 



Integral (2) represents the solution of the equation of 
radiative transfer for the case in which elastic scattering 
is described by a Gaussian phase function. 

15 The probability distribution function for the case of 

isotropic scattering, generally referred to as the diffu- 
sive limit, does not follow immediately from Eq. (2) , which 
was derived in the approximation of small angle scattering. 
However the diffusive limit can be calculated by eliminat- 

2 0 ing the angular dependence of the phase function. Thus, in 
this case P(r„, r B , T) is a path integral over the function 
Jtr (t) ] : 

By writing the trajectory in the form of a Fourier 
sine series with a fundamental period of T , an approximate 
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P(r a ,r B ,T) =f*DrU)J[z(t)} . 



expression for Eq. (3) can be derived: 
0, R>T 

PiR - T) =\(F(T)(l-^)^\ R<T < 5 > 
T 2 



T 3 B(3/2,3u s r/4-l) 



and B(a,B)is 



the beta function. This integral takes on different values 
for T<R and T>R. Note that for T<R the probability (3) is 
equal to zero, which simply reflects the requirement of 
5 causality there is not enough time for a photon initial- 
ly at A to get to point B. When R<<T, we immediately 
obtain 

This formula is a standard result of the time depen- 
dent diffusion approximation. In contrast Eq. (4) has the 
LO correct properties for times R~T. 

A much more important and interesting case for which 
the scattering is not isotropic and the phase function is 
highly peaked in the forward direction. This is the case 
of relevance to biological tissue. The solution to the 
L5 transport equation for a sample of infinite extent is giver, 
by the integral Eq. (2). Thus: 

7 

P(r A ,r B ,T)=[*~d«> [Di(t)Qexp{-[ [- — A -f ( t) 2 -io> (f ( t) 2 -i } ] dz 

J— J j 2u £ (l-g) 

(7) 

where Q is a normalization constant. The same approach 
employed in the diffusive limit can be used to evaluate 
this integral. Calculations show that the probability 
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depends upon the two factors, one resembling the "diffu- 
sion" term which dominates in the long time limit, and the 
second, which comes into play at short times, and is most 
important for the early photons that follow short, "almost 
5 straight" trajectories. The application of path integral 
techniques to the problem of light propagation has been set 
forth by Perelman et al. in "Photon Migration in Turbid 
Media Using Path Integrals", Physical Review Letters, Vol. 
72, No. 9, (1994) the contents of which is incorporated 

10 herein by reference. The application of path integral 
methods for representing light propagation in biological 
tissue is set forth in Perelman et al . , "Time Dependent 
Photon Migration Using Path Integrals", Physical Review 
(1995) , the contents of which is incorporated herein by 

15 reference. In particular, the representation for small 
angle scattering is of critical importance in accurately 
imaging embedded lesions using early arriving photons. 

To further improve upon the solution to the inverse 
problem in connection with the use of early arriving scat- 

2 0 tered photons for medical imaging. One method involves the 

calculation of the width of the most probable path repre- 
sented in Eq. (2) to define a "probe." Calculating how 
many trajectories are eliminated by blocking different 
parts of this "probe" by an absorbing object inside the 
25 tissue, one can determine the change in detected signal as 
a function of size and position of the object. Thus, the 
width of the distribution (especially for early times) and 
amount of signal obtained at the detector becomes impor- 
tant . 

30 For example, using Eq. (2) we calculated a maximum 

width of the path distribution around a classical path 
(defined as the width where path density decreases e times 
from its maximum value) for the 3 cm slab of breast tissue- 
like medium without absorption (^' s = 10 cm" 1 , and n = 0) . 

3 5 As seen from Fig. 13A, a significant improvement in resolu- 
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tion is obtained by moving the time-gate to shorter times 
(around 400 ps) . The signal level (Fig. 13B) for these 
times is around 10% of its maximum "diffusive" value, and 
for a simple shadowgram we can obtain spatial resolution of 

5 the order of 3-4 mm, which is 3-4 times better than that 
obtained using diffusive photons. 

By obtaining the distribution of paths around the 
classical path in this beam (where one estimates a uniform 
distribution with a width calculated using the path inte- 
10 gral technique, or alternatively, using a Gaussian distri- 
bution, which directly comes from Eq. (2)}, one can esti- 
mate the change in signal because part of the beam is 
blocked. This provides a relatively simple procedure for 
creating a data kernel matrix G for early arriving photons. 

5 The accuracy of this approximation and the resolution can 
be calculated from time-resolved transmission through a 
turbid medium with a screen located in the middle, as a 
function of source-detector pair position scanner across 
the medium. Fig. 14 illustrates results of calculations 

0 for the direct problem, which was performed using a beam 
with uniform distribution, and its comparison with a time 
dependent experiment. Using this simple uniform model for 
distribution of trajectories inside the beam of early 
arriving photons, one can locate the edge of the screen in 

5 a manner similar to the way it is located in a transparent 
medium. 

Using this method for more complicated geometries and 
shapes of imbedded objects, a linear system of equations 
(or G matrix representing the system) defines the temporal 
0 signal at the position of the probe as a function of the 
characteristics and positions of absorbing objects inside 
the tissue. By determining the generalized inverse G" d , 
one can establish the distribution of the absorbers inside 
the tissue. 
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In general terms the data obtained from such a 
measurement, the data vector f (t) is related to the matrix 
G F by a vector of model parameters M: f (t) = G F (t)M. The 
matrix G F (t) is provided by the operator 

G F (r,r F , t) = 

Jc(r-r F , t-t")exp(- t "~ t ' |g (r f -r„ t'-t)Q s {r g , t)dr s dt>dt" 

5 (8) 
where G(r,t) is a photon migration Green's function and Q s 
represents the spatial and temporal distribution of the 
incident laser beam. As this problem is linear, so-called 
"Maximum Likelihood" methods can be used. Details regard- 

10 ing the use of such methods can be found in Geophysical 
Data Analysis: Discrete Inverse Theory by William Menke , 
Academic Press C 1989, the contents of which are incorpo- 
rated herein by reference. 

Because the method is based on the concept of the beam 

15 of the trajectories, it provides a representation of the 
photon migration process for early arriving photons and 
thus provides the data kernel matrix in the way similar to 
methods used in x-ray tomography. The same procedure can 
be used for fluorescence measurements of tissue both in 

20 vitro and in vivo . Note that the term G(r-r F , t-t") in 

equation (8) represents the incident elastically scattered 

light, the term exp ( - t "~ t> ) represents the fluorescence, 

where r is the fluorescence lifetime and G(r F - r s , t' -t) 
represents the returning scattered light. 
25 The inverse problem can also be addressed for applica- 

tions in which no linear solution of the inverse problem 
exists which also provides the resolution and signal 
strength necessary. 
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Consider cases where a continuous distribution of 
absorbing and scattering properties inside the tissue. 
Note that absorption reduces the probability for realiza- 
tion of the particular trajectory. This can be incorporat- 
5 ed into Eq. (2) by multiplying the probability for the 
specific trajectory by a weight equal to exp ( -J> a (r) dr) , 
where the integral is calculated along the trajectory. 
Thus, Eq. (2) can be re -written in the form 



10 where the function under the integral in the exponent is 
the effective Lagrangian in a non-unif ormly absorbing 
medium. It has been generalized by introducing the func- 
tion V(r) =2/x„(r)/i a (l-g) , associated with absorption. This 
function has the same properties as a potential, and there - 

15 fore can be called the effective potential. 

Because the most probable path is disturbed by intro- 
ducing this effective potential, the problem is similar to 
finding the potential from the form of the particle's 
trajectory. This problem can be solved by using one of the 

20 minimization methods for least square problems described, 
for example, in more detail in "Generalized non- linear 
inverse problems solved using least squares criterion, " by 
A. Tarantola and B. Valette in Rev : Geophvs . Space Phys 
20, 219-232 (1982), the contents of which are incorporated 

25 herein by reference. 

In cases of non-homogeneous scattering coefficient 
H s {r) , the term in front of the action in the exponent will 
go under the integral. This term, however, does not resem- 
ble a potential, but rather a change in particle "effective 

3 0 mass" . 



P(r t ,r, 




I>r(t)expf. 




(9) 
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Retuming to systems and methods for practicing the 
invention, as described in connection with Figure 1, the 
system can be employed for measurements using -150 fs 
excitation pulses generated by a mode-locked Ti:sapphire 

5 laser pumped by a multiline argon ion laser, and a streak 
camera detection system. In another preferred method the 
incident wavelength was 800 nm, the repetition rate 76 MHz, 
and the average power 1.5 W. A small portion of the exci- 
tation light, deflected by a quartz plate to a fast photod- 

0 iode, was used as the optical triggering signal. Transmis- 
sion signals were collected by a 200 (im core diameter 
optical fiber. The other end of the fiber was imaged onto 
the streak camera slit. The fiber optic elements can be 
mounted on a rotating C-arm to position delivery and col- 

5 lection fibers as desired. This actuating mechanism can be 
operated by suitable motor and computer controlled systems 
for automatic scanning and collection. The system resolu- 
tion, 10 ps, was determined by the intrinsic response of 
the streak camera, the temporal dispersion through the 

0 optical fibers, and optical trigger jitter. 

In the case of biological tissue in which scattering 
is not isotropic and is peaked in the forward direction, 
the solution to the transport equation is given by equation 
(2) . We obtain: 



P(r f , r s , T) = 

5 where, F(T) is a normalization function. If we consider a 
slab of thickness L and introduce a coordinate system in 
which x and y axes are parallel to the slab's surface and z 
is normal to this surface, then we can approximately write, 
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P(r t ,r.,W,T) =0[ 3M3 ^ f ' g) ] T->'il-(L>+D>)/T>r"-° >T/ < x 

^ 7f J , (ID 



where D= (D x , D y ) is a vector that represents displacement 
of the detected photon from the slab's axes and the vector 
W describes the angle at which the photon emerges from the 
slab (See Figure 13A and 13B) . The components of 
5 W=(W X , W y ) are the angles it makes with the x-z and y-z 
planes respectively, and Q is a normalization constant. 

As can be seen, the probability depends upon two 
factors, one resembling the "diffusion" term, which domi- 
nates in the long time limit, and the second, which comes 
10 into play at short times, and is most important for the 

early photons that follow short, "almost straight" trajec- 
tories . 

Collected data were compared with the theoretical 
representation that follows from Eq. (11) . As shown in 

15 Figure 12, the results indicate agreement between theory 
and experiment that can be obtained by varying only one 
parameter, the scattering coefficient ^ £ . This agreement 
exists not only for the decreasing ("diffusive"! part of 
the curve but also for the initial part of the curve that 

20 represents the earliest, highly "non-diffusive" photons. 

Moreover, Figure 12 suggests that Eq . (11) has the correct 
angular dependence . 

Figure 12 shows time-resolved signal measurements for 
transmission through a slab of turbid media (Intralipid 

25 with concentration C Intr +10 ml) of thickness L=5 . 5 cm. The 
probe was displaced laterally 2 cm from the incident laser 
beam. Two curves represent two different angles of the 
probe, 6 = 4.0° and 0 = -4O° measured with respect to the direc- 
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tion of the incident beam. The measured data points are 
compared with theory (solid lines) with ji B — 3 cm" 1 , g+0.8. 

Shown in Figures 15A-15E are schematic illustrations 
of systems used to deliver and collect radiation in accor- 
5 dance with preferred embodiments of the invention. 

Figure 15A employs a single ring 5 0 in which individ- 
ual fibers 52 are aligned with an object to be imaged 60 
which can be an excised tissue sample or a portion of the 
human body in vivo. The ring can be rotated 54 to reposi- 

10 tion the fibers relative to the target and/or the entire 

ring can be scanned in a direction 56 that is orthogonal to 
the plane in which the ring 50 extends. This linear scan 
56 and rotational scan 54 of ring 50 permits the use of a 
single ring to collect high resolution three dimensional 

15 images of a lesion within the object 60. A flexible opti- 
cal fiber cable 62 connects the ring 50 to collection 
optics 64 described earlier. 

The radiation source can be coupled to the target from 
different directions. One option is to align the source 

20 along an axis that is orthogonal to the plane of the ring 
and which extends through the target 60. This can be seen 
in another preferred embodiment 70 illustrated in Figure 
15B. 

In this embodiment, a partial ring 76 is used to hold 
25 and align the fibers 52. This ring or "C-arm" 76 an be 

scanned linearly 74, rotated 75 in the plane of the C-arm, 
or alternatively, can be rotated about an axis extending 
through both ends of the C-arm. The source 72 is posi- 
tioned along an axis 77 that is orthogonal to the plane of 
30 the C-arm. The source 72 can be rigidly aligned with the 
collection array, or alternatively, can be separately 
positionable reltive to the object or the collection array. 

Separate actuators can be used to control rotation and 
translation of both the collection array and the source. 
3 5 The actuators are connected to a controller and a computer 
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programmed to control movement and position of each compo- 
nent relative to the object 60 and each other. Precise 
control of position and recording the position of each 
component is critical for image reconstruction. The com- 

5 puter uses the position data during image reconstruction. 

Another preferred embodiment is illustrated in Figure 
15C in which three rings 82, 84 and 86 are stacked to form 
a cylinder 80. More rings can also be used depending upon 
the depth of the target to be image. The entire cylinder 

0 80, or each of the individual rings 82, 84, 86, can be 
rotated 85 to provide the desired angular resolution. 

Another embodiment of each of the above systems in- 
volves incorporating a much larger number of fibers into 
each ring to provide the desired angular resolution and 

5 consequently minimize or eliminate the need for rotation. 
In the case of breast imaging, for example, in which the 
interior ring diameter is about 20 cm, and each fiber has a 
diameter of about 0.5-3 mm (preferably in the range of 1-2 
mm) , between 50 and 200 fibers can be positioned in a 

0 single ring. Relatively large diameter fibers are needed 
to provide a sufficient collection efficiency at each 
collection point due to the small number of photons being 
collected . 

In the embodiment 90 of Figure 15D, a collection 
5 geometry extending over a relatively small arc 98 is used 
along with a scanning source which can be mechanically 
moved from a first position 92 along any selected arc 
relative to the object 94 to be imaged to a second position 
96. The scanning can also be done electromechanically by 
0 coupling the laser source to any of a plurality of fibers 
positioned relative to the object 94. Such coupling sys- 
tems have been described in the above referenced incorpo- 
rated applications. Thus, a scanning optical system can be 
used to reconstruct a three dimensional image. 
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Another preferred embodiment 100 is illustrated in 
connection with Figure 15E. An endoscope or fiber optic 
probe 102 is used to deliver radiation through a first 
fiber or array of fibers 106. The distal end of the probe 

5 can be separated from the tissue 108 or can be directly in 
contact with the tissue surface. Radiation 120 from array 
106 is directed into the tissue and is scattered along 
shorter 114 or longer 112 paths through the tissue 108 to a 
position where it can be collected 114 by fiber or fiber 

0 array 104. Note that in one embodiment each fiber in array 
106 or array 104 can be used for separate delivery or col- 
lection, respectively. The retainer 110 can also be rotat- 
ed within outer housing 102 of the probe to provide differ- 
ent views of the tissue and thereby provide three dimen- 

5 sional information for image reconstruction of the embedded 
lesion of interest. 

Equivalents 

While the invention has been described in connection 
with specific methods and apparatus, it is to be understood 
0 that the description is by way of example and not as a 
limitation to the scope of the invention as set forth in 
the claims. 
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CLAIMS 

We claim, 

1. A system for detecting material within tissue 
comprising : 

5 a light source that irradiates tissue with opti- 

cal radiation along an optical path; 

an optical system that delivers the radiation 
along the optical path and that collects scattered 
optical radiation from the material; 
10 a detector that detects an intensity distribution 

of the collected optical radiation in a time interval 
between 0 and 1500 picoseconds after irradiation of 
the tissue; and 

a computer that compares the measured intensity 
15 distribution with a representation of light trajecto- 

ries in the tissue stored in a memory of the computer 
and that forms a multi -dimensional representation of 
the material within the tissue. 

2. The system of Claim 1 further comprising providing a 
20 filter to remove a portion of the collected radiation. 

3 . The system of Claim 1 further comprising an actuator 
that controls movement of the fiber optic device 
relative to the tissue and the fiber optic device 
comprises a plurality of optical fibers that collects 

25 radiation from the tissue at different non-coplanar 

angles . 

4. The system of Claim 1 wherein the detector comprises a 
with a charge coupled device. 
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5. The system of Claim 1 further comprising gating means 
for time gating the collected radiation. 

6 . The system of Claim 1 wherein the representation 

comprises a non-diffusion representation stored in the 
5 memory . 

7. The system of Claim 1 further comprising a scanning 

system to provide relative movement between the opti- 
cal system or the light source and the tissue. 

8 . The system of Claim 1 further comprising an actuator 
10 system that moves the optical system or the light 

source relative to the tissue and that generates 
position data to be recorded in the memory. 

9. The system of Claim l wherein the optical system 
comprises a ring of optical fibers extending radially 

15 from a position within the ring, the optical fibers 

having a diameter in the range between 0.5 and 3 mm. 

10. The system of Claim 1 wherein the optical system 
comprises a fiber optic probe that can be inserted 
into a body lumen. 

20 11. The system of Claim 1 wherein the collection system 

collects fluorescent or Raman scattered light from the 
tissue . 

12. The system of Claim 1 wherein the optical system com- 
prises an array of collection fibers positioned by a 
25 plurality of rings in a stacked position. 



The system of Claim 9 further comprising an actuator 
to rotate the ring in the plane of the ring. 
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14. The system of Claim 9 wherein the ring comprises a C- 
arm. 

15. The system of Claim l wherein the light source com- 
prises a laser and further comprising an optical cou- 

5 pier to couple the laser to any selected optical 

fibers of the optical system. 

16. A method of imaging material within an object 
comprising : 

irradiating the object with optical radiation 
10 along an optical path, the object having material 

located underneath a surface of the object; 

collecting scattered optical radiation from the 
material ; 

measuring an intensity distribution of the col- 
15 lected optical radiation in a time interval between 0 

and 1500 picoseconds after irradiation of the object; 
and 

forming an image of the material within the 
object from of the measured intensity distribution. 

20 17. The method of Claim 16 comparing the measured intensi- 
ty distribution with a non-diffusion representation of 
light trajectories in the tissue; and 

18. The method of Claim 16 further comprising collecting 
scattered fluorescent light from the material. 



25 19 . 



The method of Claim 16 further comprising irradiating 
the object with substantially monochromatic light from 
a laser. 
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20. The method of Claim 16 further comprising generating a 
frequency domain representation of the measured inten- 
sity distribution. 

21. The method of Claim 16 further comprising identifying 
a fluorphore within the material and determining a 
distribution of the fluorophore within the material. 

22. The method of Claim 16 further comprising detecting 
Raman shifted radiation from the material. 

23. The method of Claim 16 further comprising comparing 
the measured intensity distribution collected within 
the time period between 0 and 500 picoseconds after 
irradiation of the tissue. 

24 . The method of Claim 16 further comprising comparing a 
rise time of collected scattered fluorescent light 
with a reference, the reference being formed with a 
path integral representation of light paths through 
the object. 

25. The method of Claim 18 further comprising sensing the 
fluorescent radiation with a charge coupled device. 

26. The method of Claim 18 further comprising providing a 
plurality of optical fibers to collect the scattered 
fluorescent radiation. 

27. The method of Claim 18 further comprising moving a 
fiber optic device to deliver or collect radiation 
being directed through the tissue at different non- 
coplanar angles or positions. 
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28. The method of Claim 16 wherein the image forming step 
comprises forming a three dimensional image. 
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